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ABSTRACT 

This paper investigates a problem on the steady-state, conduct ion- convection 
heat transfer process in cylindrical porous heat exchangers. The governing partial 
differential equations for the system are obtained using the energy conservation law. 

Solution of these equations and the concept of enthalpy lead to a new approach to 
prove a theorem that the sum of inverse squares of all the positive roots of the zero 
order Bessel function of the first kind equals to one-forth. As a corollary, it is shown 
that the sum of one over pth power ( p > 2) of the roots converges to some constant. 

NOMENCLATURE 

W specific mass flow rate , Kg/m 2 .s 

C p specific heat of fluid at constant pressure , J/Kg °C 

h convective heat transfer coefficient , W/m 2 .°C 

R radius of cylindrical board , m 

r radial coordinate , m 

z axial coordinate , m 

T s solid temperature , °C 

T s o solid temperature at circumference , °C 

Tf fluid temperature , °C 

Tf o inlet fluid temperature , °C 

A WCpR/Xr , dimensionless parameter 

B / 3hR/WCp , dimensionless parameter 

‘ (T, - T„)/(T m - Tj 0 ) , dimensionless temperature of solid media 
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T (T, - T, a )/(T. 0 - T/o) , dimensionless temperature of fluid 
x z /R , dimensionless axial coordinate 
y r/R , dimensionless radial coordinate 

Greek Letters 

/? specific area of heat transfer , m 2 /m 3 
7 A,/A r , ratio of the axial and radial thermal conductivity 
X z axial thermal conductivity of solid , Wjm.°C 
A r radial thermal conductivity of solid , W/m°C 


INTRODUCTION 

Porous heat exchangers play a significant role in many engineering appli- 
cations, such as cryogenics, thermal storage systems, and chemical reactors. Such 
systems lead to a set of Partial Differential Equations (PDEs) with a strong coupling 
between equations for the solid and the fluid phases. Analytical schemes to solve 
a general class of problems in this area include the method of separation of vari- 
ables (refs. 1 and 2), Riemann method (ref. 3), orthogonal collocation techniques 
(ref. 4) and collocation-perturbation scheme for packed beds (ref. 5). Siegwarth 
and Radebough (ref. 6) present a numerical technique to solve the above problems 
with variable physical properties. Lin, Guo, and Wang (ref. 7) present a com- 
bined orthogonal collocation-perturbation method to solve the temperature field in 
a cylindrical porous heat exchanger. 

In this paper, the physical problem presented in (Ref. 7) is reconsidered. 
Using the approach of separation of variables, the PDEs associated with the prob- 
lem are reduced to two Boundary Value Problems (BVPs) of Ordinary Differential 
Equations (ODEs), Solution of the resulting BVPs leads to a new method to prove 
a theorem regarding roots of the zero-order Bessel function of the first kind. As 
a corollary, it is shown that the sum of one over pth power (p > 2) of the roots 
converges to some constant. 

The theory of Bessel functions and related functions is well established and 
the main results are summarized in textbooks and mathematical manuals (refs. 8 
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to 10). The properties of Bessel functions concerning operations of differentiation 
and integration with respect to their order have also been studied by Apelblat and 
Kravitsky (ref. 11). However, little is known about the convergence and summation 
of series for roots of Bessel functions due to no explicit expressions for the roots 
except for ±1/2 order Bessel functions. Although in 1874 Rayleigh derived a gen- 
eral form of the theorem discussed here by applying an Euler’s formula (ref. 8), 
his approach is suitable for even power series only. Furthermore, he did not study 
the convergence of such series for a general case. Therefore, research is needed to 
further understand the properties of Bessel functions. The new method to prove the 
theorem stated above is based on the mathematical model of a cylindrical porous 
heat exchanger. A brief description of this model is presented next. 

MATHEMATICAL MODELING 

The model considered here is a semi-infinite cylindrical porous heat ex- 
changer as shown in (fig. 1). Let r and z be the radial and the axial coordinates, 
and R be the radius of the cylinder. A fluid flows through the porous media in axial 
direction from left to right. The inlet temperature of the fluid is T f o which can be 
higher or lower than the board circumference temperature T, 0 . Let W, C p , and h 
be the specific mass flow rate, the specific heat of fluid at constant pressure, and 
the convective heat transfer coefficient respectively. 

Porous materials used in such applications exhibit anisotropic behavior. It 
is assumed that the thermal conductivity of the solid is symmetric about the axis 
of the cylinder. Let be the specific area of heat transfer , k z and k r be the axial 
and the radial thermal conductivities of the solid, and 7 be the ratio of k z and k T . 

In the derivation to follow, we assume that: 1) the physical properties and 
convective heat transfer coefficient between porous substance and fluid are con- 
stants, 2) the dimensions of the porosity and solid particles are very small com- 
pared to the overall dimension of the heat exchanger, and therefore, the porous 
material can be treated as continuous media, 3) the fluid thermal conductivities 
are negligible compared to the solid thermal conductivities, and 4) the inner wall 
temperature is kept constant, and there are no thermal resistances between wall 
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and porous media or fluid. Finally, for simplicity, it is considered that T . o > Tf 0 . 
However, the formulation presented here equally holds for T* < T/o. 

The differential equations of the system may be obtained by applying the 
energy conservation law to a micro unit volume (ref. 7). Nondimensional forms of 
these equations and the boundary conditions are given as follows: 

Solid Phase : 

d 2 t ldt_ _ dT 

dy 2 + y dy ^ dx 2 dx ^ 

Fluid Phase : 

f = B( t -T) (2) 

Boundary Conditions : 


[* = i, y = 1 
< r = o, x = o 

t = 1, x — + +oo 


(за) 

(зб) 
(3c) 


where f(= (T, — T/ 0 )/(T J o — T/ 0 )) and T(= (Tj —Tf 0 )/(T l0 — T/ 0 )) are dimensionless 
temperatures of the solid and the fluid, x(= z/R ) and y(= r/R ) are dimension- 
less radial and axial coordinates, and /!(= WC p R/k r ) and B(= 0hR/WC p ) are 
dimensionless parameters. Here T, and Tj are the solid and the fluid temperature 
distributions. It should be noted that an additional boundary condition is required 
to completely define the problem. However, it is not needed here and therefore it 
is not considered. 

Equations (1) to (3) can be used to find the temperature distributions of the 
solid and the fluid phases. It will be shown that these equations can also be used 
to prove the following theorem: 

Theorem: Let a n denote the nth positive root of the zero-order Bessel function 

of the first kind, then 



2 2 


( 4 ) 


Corollary: For p > 2, series converges to a constant C(p ) which 

depends on p. 
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Proof: In order to prove the above theorem, we begin by eliminating T from Eqs. 
(1) and (2). This leads to 


Let 


d d 2 t IcH 
dx dy 2 + y dy + 7 


dx 2 




3H 

dy 2 


1 dt d 2 t 
+ yd^ + J d^> 


( 5 ) 


t=l-XY (6) 

where X and Y are functions of x and y respectively. Using Eqs. (3) and (5) 
and the method of separation of variables, we obtain the following two equivalent 
boundary value problems of ordinary differential equation: 


r y" + l Y '-\y = o 

' Ym 0, y = 1 
k \Y\ < oo, y = 0 


and 


(7a) 

(76) 

(7c) 


7 X‘" + ~{BX" + (A - AB)X' + A BX = 0 (8a) 

x = 0, x ► +oo (86) 

where a prime( / ) on X (Y) represents the derivative with respect to x (y), and 
A is a separation constant. Equation (7c) suggests that Y is bounded at its center. 
Equation (7) can be brought to standard Bessel equation form by a simple linear 
transformation. This is a general Sturm - Liouville eigenvalue problem (ref. 14). 

It has nontrivial solutions only when A n = — a^(a n ^ 0). Using the properties of 
Bessel functions, the solution of Eq. (7) for any ot n may be written as 


Y n = c n J 0 (a n y ) ( 9 ) 

where Jo(x) is the zero order Bessel function of the first kind, a n is the nth positive 
root of /„(*), ^d c n is an undetermined constant. 

Let the trial solution of Eq. (8) be given as 

X = a'e* (10) 
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Substitution of Eq. (10) into Eq. (8) leads to the following characteristic equation, 

7/? 3 + ~fB0 2 - (a 2 + AB)j3 - a 2 n B = 0 (11) 

Equation (11) is cubic in j3. The discriminant A of this equation is given as 
A = — ^-r{4a 3 (B 2 7 - a 2 ) 2 + 18A£? 2 a 2 7 + AB 3 (AB 4- 2a 2 )+ 
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4AB(A 2 £ 2 + 3a* +3a 2 AB)} (12) 

Equation (12) suggests that A < 0. Thus, using the theory of cubic equations, it 
follows that all roots of Eq. (11) are real (ref. 12). Let these roots be given as 9ni, 
P n2 , and 0n3- From Eq. (11), it follows that 

9nl + 9n2 + 9n3 — < 0 ( 13 ) 


and 

9nl ' Pn2 • PnZ = Q n-®/7 > 0 (1^) 

Equation (13) suggests that atleast one of the roots is negative and Eq. (14) indi- 
cates that negative roots appear in pair. Thus, it is concluded that Eq. (11) has 
two negative roots and one positive root for each a n . From physical consideration 
(or Eq. (8b)), the positive root is disregarded. The general solution for X may now 
be written as 

X, = £ a n.' e/3mX (15) 

1=1 

Substituting Eqs. (9) and (15) into Eq. (6), we get 


t = 1 - E M a nV) ■J2a ni e 0n,x 


71 = 1 


»= l 


where a ni = c n a' ni . From Eqs. (1), (2) and (16), we obtain 

2 

i 

n=i t=i 

Differentiation of Eq. (17) with respect to x yields, 


1 00 

AB 


T = 1 — rrr ■ E Jo(a n y ) • E a m^ n,x ■ (*» - 7^n, + AB) 


dT 

dx 


AB 


£ MW ) ■ £ a niPni • («n " 7 Pli + AB)e 


t Pm X 


(16) 


(17) 


(18) 


71=1 


*=T 
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Using Eq. (17), the inlet condition (Eq. (3b)), and the property of orthogonality 
of J 0 (a n y ) in the interval [0, 1], we derive 


£ ani(o£ - 7 Pli + AB) = ■ 2AB 

j'=l &n * 

where J\(x ) is a Bessel function of the first kind of order one. 


(19) 


ENTHALPY RELATIONS 

Some enthalpy relations help prove the above theorem. Let H(x ) and H(x -f- 
dx ) be quantities of enthalpy carried by the fluid across the axial planes at x and 
x + dx respectively (Figure 1). Expression for H(x) is given as 


H(x) = WC P • f R {rdr • /” d9 • [T /0 + T(T, 0 - T /0 )]} 

JQ J 0 


( 20 ) 


A closed form expression for H(x) can be obtained by substituting Eq. (17) into 
Eq. (20), and integrating the resulting equation. Expanding H (x +dx) in a Taylor’s 
series, and neglecting the second and the higher order terms, we obtain 


dH(z) = H(x + dx) - H(x) = ■ dx 


( 21 ) 


where dH(x) is the enthalpy variation of the fluid passing through the control 
volume between x and x + dx (Figure 1). The total change in enthalpy of the fluid 
is obtained by integrating Eq. (21) over the length of the heat exchanger. Thus, 

A Hr-l + ~dH(x) = l*~ d JM. dl (22) 

Using Eqs. (18), (20), and (22), we obtain 

_ ^C f T, A ±{J ^ _ ^ + AB)] (23) 

n=l 


t=l 


AHq°° can also be evaluated directly by subtracting the enthalpy of the fluid at 
the inlet (= H{ 0)) from that at the outlet (= H( oo)). Expressions for H( 0) and 
H( oo) Eire given as 

H( 0) = WC p kR 2 -T /o 

and 

H{ oo) = WC p xR 2 ■ Zo 
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Hence, 


A £T 0 +o ° = WC p irR 2 ■ (T, 0 - T /0 ) 


Comparing Eqs. (23) and (24), we obtain 

— + = l (25) 

AU n=l i=l 

Finadly, substituting Eq. (19) into Eq. (25) and simplifying, we obtain 

+“ i 1 

V — = - 

A* a 2 2 2 

n=l a n * 

This proves the theorem. Since function Jo(y) is symmetrical, it follows that 

+o° 1 I 

n=— co n * 

In order to prove the corollary, observe that all positive roots of Jo(x) are greater 
than 1 (ref. 8). This implies that for p > 2 

-V < 4 ( 27 > 

an al 

Equations (4) and (27) suggest that the series 

°n(p) = E ^ 

t=l 

increases monotonically and it is bounded. Thus, by monotone convergence theo- 
rem (ref. 13), it follows that the series C n (p) (n — ► oo) is convergent. 


CONCLUSION 

A mathematical model for a steady-state conduction-convection heat transfer 
processes in a cylindrical porous heat exchanger has been investigated. It has been 
shown that the equations resulting in this model may be used to prove a theorem 
and a corollary regarding roots of the zero-order Bessel function of t he first kind. 

Research reported in this paper provides an alternate approach to prove a 
theorem in this area. The advantages of this physical approach are that if enriches 
physical understanding of a theorem, and that we may avoid difficulties emerging 
from rigid mathematical arguments. 
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Figure 1 . Porous heat exchanger model 


498 



